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^ ■ Abstract 
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■rj" ' tum information processor are presented, and shown to be consistent with a detailed model of the 

system-plus-apparatus used for the experiments. The quantum operation studied was the quantum 
Fourier transform, which is important in several quantum algorithms and poses a rigorous test for 
the precision of our recently-developed strongly modulating control fields. The results were ana- 
C I lyzed in an attempt to decompose the implementation errors into coherent (overall systematic), 

^ I incoherent (microscopically deterministic), and decoherent (microscopically random) components. 

^^' This analysis yielded a superoperator consisting of a unitary part that was strongly correlated with 

the theoretically expected unitary superoperator of the quantum Fourier transform, an overall at- 
^ ■ tenuation consistent with decoherence, and a residual portion that was not completely positive - 

although complete positivity is required for any quantum operation. By comparison with the results 
of computer simulations, the lack of complete positivity was shown to be largely a consequence 
of the incoherent errors during the quantum process tomography procedure. These simulations 
further showed that coherent, incoherent, and decoherent errors can often be identified by their 
distinctive effects on the spectrum of the overall superoperator. The gate fidelity of the experimen- 
tally determined superoperator was 0.64, while the correlation coefficient between experimentally 
determined superoperator and the simulated superoperator was 0.79; most of the discrepancies 
with the simulations could be explained by the cummulative effect of small errors in the single 
qubit gates. 
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I. INTRODUCTION 



In order to develop larger and more powerful quantum information processing devices, it 



is essential to quantify the precision with which they can be controlled. 



'his information 



is generally reported as a single number, the fidelity of the operation jl|, |2, |3|. Although 
fidelity is a reasonable measure of control, it gives experimentalists little useful informa- 
tion about what went wrong or how to irnprove their control over the quantum system. 
Quantum process tomography (QPT) p, ly, |7| provides additional information that may be 
useful in this regard, by yielding an estimate of the quantum operation that was actually 
implemented. This in turn provides a stringent check on the completeness and accuracy of 
the system-plus-apparatus model used to design the implementation. It is nevertheless a 
challenging task to interpret the deviations of this estimate from the superoperator implied 
by the model in terms of specific defects in the model. Additional errors introduced during 
the QPT procedure itself further complicate the analysis of the QPT results. In this paper 
we explore these issues via a concrete case study, in which QPT is performed on a previ- 
ously reported three-qubit quantum Fourier transform (QFT_) implemented using a nuclear 
magnetic resonance (NMR) quantum information processor 8|. 

The dynamics of an isolated quantum system are described by the Schrodinger equation, 
which gives rise to an A^ x A^ unitary operator, where A^ is the dimension of system's Hilbert 
space. Open quantum systems, however, generally interact with an inaccessible environment 
and thereby undergo decoherence p|. Furthermore, in the case of expectation value mea- 
surements as in liquid-state NMR, each element of the statistical ensemble may undergo a 
slightly different unitary operation which, though nonrandom, is difficult to distinguish from 
decoherence J9|,|l0|. The statistics of measurements on open quantum systems are generally 
described by an A^ x A^ density operator, and the evolution of a density operator under an 
incoherent distribution of Hamiltonians and/or interactions with an environment, although 
non-unitary, remains linear and is described by a superoperator. The goal of quantum 
process tomography is to determine this superoperator. 

Methods for implementing QPT have been presented Refs. p. If , m Il2l |l3[, and two- 
qubit NMR implementations of QPT have previously been reported IJ, |l5| . In the present 
paper, QPT is carried out with the aim of validating the mathematical model of the system- 
plus-apparatus used to design the NMR implementation of a multi-qubit, entangling unitary 



operation, and to identify the types and strengths of the errors that occurred. We have found 
it useful to classify the errors as coherent, incoherent, and decoherent, because each class is 
related to specific short-comings in the experimental implementation. Coherent errors are 
systematic errors in the net unitary operation that was actually implemented. Incoherent 
errors refer to unwanted unitary evolution which is not uniform across the ensemble of spin 
systems in the NMR sample, so that even though each member of the ensemble undergoes 
strictly unitary evolution, the evolution of the ensemble averages appears non-unitary. The 
effects of incoherence are reversible, at least in principle, and knowledge of the coherent and 
incoherent errors can be used to design better quantum gates (unitary operations). In con- 
trast, decoherent errors are due to unknown interactions with an inaccessible environment, 
so they are not reversible and can be eliminated only by relatively costly changes to the 
apparatus or the way in which information is encoded within it. 

An important benefit of QPT is that it poses a rigorous test of the accuracy of the math- 
ematical model of the system-plus-apparatus used to design and interpret the experiments. 
This is done by comparing the experimental results of QPT to the results of computer sim- 
ulations of the complete QPT procedure, based upon this same model. Simulations based 
on the model used here were also used in designing the strongly modulating control fields by 
which both the desired unitary operation as well as all the unitary operations needed for the 
QPT procedure were implemented |lfl] . This was done by minimizing the difference between 
the desired qubit rotation operator and the quantum operation obtained by simulating the 
effect of the control fields on the spins in the molecule used for the experiments. As a result, 
any incorrectness in the model directly affects the reliability of the experiments, but in ways 
that, by definition, differ from the simulations. This can suggest ways to improve the model, 
after which further simulations will pinpoint the remaining experimental errors. 

The operation on which QPT was performed is the quantum Fourier transform (QFT). 
The QFT constitutes a key subroutine in several quantum algorithms [la, 113 as well as 
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in certain methods for simulating quantum dynamics on a quantum computer |18l . I19I |. In 
algorithms such as Shor's factoring algorithm, the QFT is used to extract periodic features 
of wave functions, while in simulations of quantum dynamics it is used to move between the 
position and momentum representations. The QFT is defined as follows: 

f/QFT|x) = ^5^e2-^-'/^|x'). (1) 



x'=0 



The QFT has been expressed J20|,|21[ as a sequence of one-qubit Hadamard gates Hj, which 
rotate the jth qubit from a computational basis state to an equal superposition of compu- 
tational basis states, and two-qubit conditional phase gates Bjk{6), which rotate the phase 
of qubit A; by ^ if qubit j is in the state |1). In this notation, the complete gate sequence of 
the three-qubit QFT is (reading from right-to-left) 

SwapnH^B^si^) B,,{^) H^B^^ifj ^i , (2) 

where Swapjk is a swap gate between qubits j and k (see Fig. [T]). This gate sequence has 
)een implemented via NMR to demonstrate the ability of the QFT to extract periodicity 



and as part of the quantum baker's map 



22| 



The remainder of this paper is organized as follows. Section II gives an overview of the 
experimental and computational procedures used for QPT, together with the metrics by 
which the results were compared to those of the theoretical model. Section II describes 
incoherent errors and their effects on QPT, in particular the apparent lack of complete posi- 
tivity of the results, while section III describes the experimental system and procedures used 
to implement QPT in detail. Section IV presents a complete description of the system-plus- 
apparatus model by which the results of QPT were interpreted. This is followed in section 
V by an overview of the experimental results, and in section VI by a detailed comparison 
of the results with the model's predictions. Finally, section VII contains an analysis of the 
discrepancies between the experimental and simulated results, with the goal of determining 
their probable origins. The paper concludes with brief discussion of the implications of our 
findings for the field of quantum information processing as a whole. 

II. QUANTUM PROCESS TOMOGRAPHY 

There are several methods of performing QPT outlined in the literature. Some of these 



methods 



B 



12| | require increasing the Hilbert space size beyond that of the system whose 
dynamics are to be studied. This is unappealing for current experimental studies of quantum 
information processing where qubits are at a premium. The procedure used in this work (see 
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Fig. 121) is similar to those given in Refs. p, [^ , and results in an A^^ x A^^ complex- valued 
matrix, hereafter referred to as a "supermatrix." Because of the inevitable experimental 
errors made in the QPT procedure, this supermatrix will not, in general, correspond to any 



completely positive and trace-preserving superoperator (ergo quanturn operation), unitary 
or otherwise. Therefore, we use methods recently introduced by Havel [isl |l5| to obtain the 
best least-squares fit to the measured supermatrix that does correspond to such a quantum 
operation. This is expected to yield a better estimate of the superoperator that was actually 
implemented, since it will subsequently be shown that the non-completely positive part of 
the superoperator is largely due to errors committed during the readout steps. 

The QPT procedure in Fig. |21 requires implementation of the desired unitary operation 
on a complete set of A^^ known input states, followed by determination of all of the resulting 
output states via quantum state tomography. Once all of the experimental input and output 
states have been completely determined, and given that the input states have been chosen 
to be linearly independent, the experimental supermatrix can be calculated as 

Mobs = 3?out 3^in ) (3) 

where %„ and CRout are supermatrices whose columns are the "columnized" experimental 

input and output density matrices col{p\„) and col{pop), respectively, as determined by state 

tomography (see below). In quantum information processing, the Hilbert space is usually a 

tensor product of the 2-dimensional Hilbert spaces of its constituent qubits, so that N = 2^ 

grows exponentially with the number of qubits n. As a result, QPT can in practice be carried 

out on only a few qubits at a time. Nevertheless, even a large quantum computer is expected 

to be based upon such local operations, so a complete analysis of small implementations is 

a prerequisite to understanding the issues involved in the control of larger systems. 

We now introduce the measures of state and gate fidelities that were used to summarize 

the results of QPT. The accuracy with which the initial states were created was quanti- 

tated by the correlation between the desired input state and the one determined by state 

tomography iQi], i.e. 

^f \ tr(pthPin) , .^ 

C(pth,Pin) = -=^==^, (4) 

Vtr(p^^Jtr(pfJ 
where pth and p\n are the traceless parts of the density matrices of the desired and measured 
input states, respectively. The input states were prepared from the equilibrium spin state 
by means of suitable non-unitary operations. Sin, and the magnitudes of their traceless parts 
taken as the reference against which all subsequent losses of coherence (or magnetization in 
NMR) due to non-unitary evolution were measured. This was done by scaling the correlation 
of the output states by a factor which measures the loss of coherence jsf, yielding the so-called 



attenuated correlation 



/tr(p2 ) 
CA(Pth,Pop) = C{pth,Pop)\h^rTV ' (5) 

where pth and pop are the traceless parts of the density matrices of the theoretical and 
measured output states, respectively, and p\„ that of the corresponding input state from 
above. 

The average state correlation and average attenuated state correlation are reasonable 
measures of the overall fidelity of the implemented quantum operation, but a measure that 
is more clearly independent of the choice of initial states is the correlation between the 
superoperator matrices, i.e. 



C(Sth, Sop) = tr(SjhSop)/y'tr(§{,§th)tr(§ip§op) . (6) 

Since this quantity is not sensitive to the overall loss of magnetization, we will also use the 
attentuated correlation between supermatrices, namely 

CA(Sth,Sop) = tr(S,VSop)/tr(§tV§th) = F,i§l§,,) , (7) 



where F^ is the entanglement fidelity defined by Schumacher jl|. Accordingly, we will refer 
to this quantity as the gate fidelity. It can be shown that this fidelity satisfies CA{pth, Pop) > 
Cyi(§th, §op), with equality if ti^pf^) = tT{p^^^) for all the input states [2|- 

Our analysis utilizes two different representations of quantum operations, the supermatrix 
representation and the Kraus operator sum representation. The supermatrix S operates on 
the density matrix as 

C0/(pfin) = S C0/(pini) , (8) 

where p\r,\ and pfjn are the initial and final density matrices, and the "col" operation stacks 
the columns of the density matrices on top of one another in left to right order. This results 
in an A^ X A^ matrix becoming a column vector of length A^^ . We will primarily work in two 
different supermatrix bases, the computational (or Zeeman) basis and the product operator 
basis. In the computational basis the rows and columns of the 2" x 2" density matrix 
are labeled by the binary expansion of their indices from |0 . . . 0) to |1 . . . 1). The product 
operator basis is related to the structure of the NMR Hamiltonian and rows and columns of 
the supermatrix in this basis are labeled as /2 ® ... ® /2, ^2 ® ••• ® -^2 ® ctx, to cr^ (g) ... ® Cz, 
where "®" denotes tensor multiplication, the cr's are the standard Pauli matrices, I2 is the 
2x2 identity matrix and each term has n factors (see Ref. 2^ for details). 



The superoperator of a quantum operation can also be expressed as a Kraus operator 
sum. This sum is comprised of a set of A^ x A^ matrices, Ak , such that 

Pfin = ^^Akp^Al (9) 

where the sum may require as many as A^^ terms. It is easily seen that a Kraus operator 
sum preserves the Hermiticity of p\r,\, and that it preserves the trace if and only if 

J2A^' = 1- (10) 

If the operation is purely unitary there is only one Kraus operator, which is just the unitary 
operator in question. In general, however, there are an infinite number of equivalent Kraus 
operator representations of a given quantum operation. The least-squares fitting procedure 
described below yields a Kraus operator sum representation in which the Kraus operators 
are both mutually orthogonal, minimum in number and sorted by the size of the contribution 
each makes to pf\„. 

III. INCOHERENT ERRORS 

A superoperator is completely positive if it admits a Kraus operator sum representation. 
Although this condition holds for the physical processes usually studied in the quantum me- 
chanics of open systems, including unitary operations and decoherence via weak interactions 
with a Markovian environment, it is not difficult to construct situations in which not only 
complete positivit y, b ut even the distinction between a state and its transformations breaks 



down 



n 



m 
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26l l27l |. Such situations can only arise when the initial state of the system 
is not pure, but can be represented by a probability distribution over an ensemble of pure 
states with density matrix p\„i. Suppose the ensemble's probability distribution depends on 
some classical parameter c (usually space or time) and that pini also represents an average 
over c (see below for a concrete example). Then if the applied transformation % likewise 
depends on c, so that the transformation is correlated with the states in the ensemble, then 
the final density matrix pfjn will generally not be equal to the result of applying the average 
of % over c to p\r,\. Indeed there will usually be no superoperator which maps every pos- 
sible Pini to the correct pfin. In principle one could define a non-linear transformation that 
produces this mapping, but this is not very useful in practice because the large amount of 



information needed to define a general non-linear transformation is neither readily available 
nor easy to work with. 

Incoherent errors are precisely those which depend upon some classical parameter that 
labels the members of an ensemble, and so generate a correlation between the states in the 
ensemble and this parameter. In the NMR experiments described below, the main source 
of incoherent errors is the spatial inhomogeneity in the RF (radio-frequency) field over the 
sample volume. This may be expressed by writing the Hamiltonian for the interaction of 
the spins with the RF field in the form 

H^,{t;f) = a(t)75i(f)e-^^-*-/V,e^^-*-/2^ (11) 

where a{t) is the envelope of the field, 7 is the gyromagnetic ratio of the spins, Bi{f) is 

the magnitude of the field at the point r , Au is the difference between the frequency of the 

RF field and the resonance frequency of the spins in the static magnetic field Bq, and the 

(t's are the usual Pauli matrices. Thus, the unitary operation U = exp ( J dt //rf) depends 

upon the position r in the sample as well as time, i.e. U = U{t;f), as does the resulting 

spin state p{t]r) = U{t]r)p{0;f)W{t]f). The equilibrium state of the spins in the static 

field Bq is independent of r over the sample volume, so we may take p(0; r ) = p(0) to be 

constant. The density operator that is observed is the integral over the sample volume V, 

i.e. 

p{t)= / dfp{t;f)= / dfU{t;f)p{0)U\t;f) = §(t)p(0). (12) 

Jv Jv 

where §(t) denotes the net superoperator of the actual quantum operation implemented. 
We can express this in matrix form as [l3| 

co/(S(t)p(0)) = (^ f drU{t; f) U{t- f )) co/(p(0)), (13) 

where U is the complex conjugate of U and the ^^coV^ operator maps N x N matrices to 
(A^^)-dimensional column vectors, as described above. 

Describing the evolution of an ensemble by integrating over its spatial degrees of freedom 
works well for a single operation. When a second operation §2 is applied following the first 
§1, however, we must take the spatial correlations produced by the first into account in 
computing the correct overall superoperator. In terms of the foregoing equations, this may 
be expressed as 

S2S1 = /" drl (1/2(^2) ® f/2(f2)) f dn (y,in) ® u.in] 

8 



= J df,J dn (f/2(r2)f/i(r-'i)) ® (f/2(r2)f/i(fi)) (14) 



where § is the true superoperator of the combined operations. Clearly, a similar result 
could be obtained if we were integrating over time, although that is not relevant to the 
experiments described in this paper. A simulation of an example of the lack of composibility 
of superoperators describing incoherent errors is shown in Fig. El 

It follows that the A^^ x A^^ supermatrix produced by our QPT procedure (Fig. ^ cannot 
be expected to precisely correspond to any physical process, and it will usually even fail to 
define a completely positive and trace-preserving superoperator. In this event, an improved 
estimate of the supermatrix of the quantum process that is realized by the QFT implemen- 
tation can be obtained by making the smallest possible change to the supermatrix derived 
from QPT so as to make it completely positive and trace-preserving (CPTP). 

One can, of course, always convert a Kraus operator sum into an equivalent supermatrix 
via the ^^col" operation introduced above. 

To go the other direction, i.e. to take a supermatrix S and convert it into an equivalent 
Kraus operator sum, we first rearrange its elements into a Hermitian supermatrix called the 
Choi supermatrix |l3l |. 

N-l 

7= ^(E,,®/jv)§(/iv®^ii), (16) 

ij=0 

where In is the N x N identity matrix and Eij is a matrix of zeros except for a 1 in the 
(i, j)th position. Then, if T = J2k ^kVkvl ^^ ^^^ spectral decomposition of T and we assume 
that all the eigenvalues Xk > 0, our Kraus operators Ak may be shown to be 

col{Ak) = ^/XkVk (fc = l,...,Ar). (17) 

It follows from the above that a superoperator S is completely positive if and only if its 
Choi supermatrix is positive semi-definite. In the event that an experimentally determined 
supermatrix Sobs does not have a positive semi-definite Choi supermatrix, it has been shown 
|l3| that the completely positive superoperator §' closest to it in the least-squares sense may 
be obtained simply by setting the negative eigenvalues of its Choi supermatrix to zero to 



get a new matrix T', and mapping it back to §' via Eqn. p6p with § and T swapped. This 
procedure, however, will result in a supermatrix §' that does not preserve the trace of the 
density matrix. That condition can be reimposed by subtracting 

AT' = N-^coI{In)coI\In)7' - Inxn (18) 

from 7', with the result that T' is no longer positive semi-definite. It has further been shown, 
however, that iterating on these two procedures generates a sequence of supermatrices which 
converges to the Choi supermatrix corresponding to the superoperator that is closest to S 
and is both trace preserving and completely positive. This is described for the experimental 
results in the Section IVIIII 

In order to quantify the extent to which a supermatrix violates the completely positive 
requirement, we define the positivity as the ratio of the sum of eigenvalues of the Choi 
matrix corresponding to the supermatrix over the sum of positive eigenvalues of the said 
Choi matrix. For a completely positive superoperator the positivity is equal to 1, while the 
presence of negative Choi eigenvalues causes the positivity to be less than 1. 

IV. THE SPIN SYSTEM AND THE EXPERIMENTS 

The exp eriments were implemented on a three-qubit NMR quantum information proces- 
sor |28Ll29J. The three qubits used were the three carbon spins in molecules of ^^C-enriched 
alanine in an aqueous solution. The internal Hamiltonian of this system has the form 

H, = vrX:^^^b + iE^c;c^c-^C' (19) 

«=1 j>i=l 

where Uq are the Larmor frequencies of the spins and Jqq are the strengths of the couplings 
between them, both in frequency units. In our indexing scheme, a superscript 1 labels the 
carbonyl carbon of alanine, 2 the alpha carbon and 3 the methyl carbon. A separate time- 
dependent external Hamiltonian, shown in Eqn. (jllll . must be added whenever an RF pulse 
is applied to rotate qubits. In the 7 Tesla magnet used for the experiments, the resonant 
frequency of carbon-13 is approximately 75.468 MHz. Frequency changes, also known as 
chemical shifts, among the spins introduce differences Vq — Vq = 9456.5 Hz, Vq — Vq = 2594.3 
Hz and ^q — t-'q = 12050.8 Hz. The coupling constants between the three spins are Jqq = 
54.2 Hz, Jq'q = 35.1 Hz, and Jq'q = — 1.2 Hz. In the absence of RF, the Ti relaxation times 
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of the three spins are all longer than 1.5 s, while the T2 relaxation times are longer than 400 
ms (see Table H] for exact numbers). 

To rotate qubits we used the strongly modulating pulses introduced in Refs. J3|, |lO[ . The 
first generation of strongly modulating pulses were designed to perform the desired propaga- 
tor while refocussing the spins' known internal Hamiltonian jSj; the second generation was 
designed to also compensate for RF field inhomogeneity [lO|. This was done by including the 
RF inhomogeneity profile in the simulations and, hence, in the target function minimized 
to design the compensated pulses. Pulses designed in this way performed near-optimally 
over the range of RF powers experienced by the ensemble of spins. While the simulated 
peak fidelities were lower than with the uncompensated pulses, the actual performance of 
the compensated pulses in the spectrometer was greatly enhanced. 



The method o 
implementations 



iraplementing the QFT via NMR is the same as in our previous QFT 
|22J . with the added benefit of the RF pulses that compensated for 
RF-inhomogeneity. The pulse sequences for the Hadamard and conditional phase gates are 
derived from an idempotent or correlation operator description of the propagators |2i| • The 
Hadamard gate pulse sequence is: 

This pulse program reads: apply a pulse that rotates spin j by 90° about the y-axis, followed 
by a pulse that rotates j by 180° about the x-axis. The Bj^ gate, can be implemented using 
the coupling between qubits and the following pulse sequence: 



■''-'4)-w:-(i):-(r-(iV' » 



y \ / X ^ ^ y 



where (f) is an arbitrary phase. The notation ( ^Jj j represents a time interval during which 
cr^cr^ evolution occurs while chemical shifts and all other qubit couplings are refocused, while 
the superscript j, k denotes a pulse which rotates only spins j and k. The final three pulses 
in the above sequence perform a rotation around the z-axis. In our previous work [8] the 
bit reversal was implemented experimentally. Here the bit reversal is achieved by simply 
renaming the bits. 

We stress that the various building blocks of the QFT implementation, i.e. the individual 
pulses and evolution periods, remain unchanged for all experiments. After creation of each 
initial state a fixed sequence is applied that is independent of the initial state. 

11 



The input states, {p\n}, used were the 64 product operator states (neglecting the large 
but undetectable identity component). These are (actually the identity Jg), and states 
such as o"x, (7^, and <Jy<J^ where the superscripts represent the spin indices, as above. These 
states were chosen because they are orthogonal and easy to create on a liquid-state NMR 
system. We order the states lexicographically starting with 0, followed by (J^ = I2®l2® <^x, 
cr^ = /2 ® CTx ® /2, and so on until crlo'lo''^ = cTz (8 0"^ ® CTz (the complete ordering may be found 
in Fig. IH}. The QFT gate sequence described above is applied to each of these 64 states to 
obtain the corresponding set of output states {pop}- 

In any given three-spin spectrum only 24 of the 64 product operators are observable. 
Thus, to completely reconstruct the state of the system, each experiment is repeated 7 
times with different readout pulses appended to the experiment. The real amplitudes of 
the peaks in the spectrum were then used as the coefficients of the product operators which 
were transformed into the corresponding observable product operators by the readout pulses. 
This procedure, known as state tomography [2|, was done for all input and output states. 

V. THE SYSTEM-PLUS- APPARATUS MODEL 

Quantum process tomography serves experimentalists in two ways. First, it provides a 
thorough test of the model of the system-plus-apparatus used to design and interpret the 
experiments. Second, it provides a complete picture of the net effects of all the errors made 
in implementing the experiment and, in conjunction with the model, allows these effects to 
be unravelled into specific defects in the apparatus - and perhaps also the model. Before 
attempting this latter task, therefore, we wish to fully present, discuss, and justify the model 
we used for QPT on the QFT. This model includes the following attributes: 

• The Hilbert space of the system in question and its Hamiltonian. 

• The Hilbert space of the larger system with which the given system interacts coher- 
ently, and the Hamiltonian describing these interactions. 

• The relaxation superoperator of the system, and some knowledge of how the larger 
system relaxes. 

• Bounds on the precision of the classical control fields applied. 
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• The distribution of incoherent variations in these fields across the sample. 

Both the fields and the positions of the spins can be treated classically. Since the larger 
system contains only a total of eight spins, we can simulate it exactly. Further information, 
for example the correlation times or full spectral densities of the noise generators driving 
relaxation J3^, could have been included in the model, but the above proved adequate to 
explain most of the experimental observations for the QFT. 

The first attribute is the Hamiltonian for the three carbon- 13 spins of alanine used as 
qubits for the QFT, which has the form of H^ given in Eqn. (jl9|) . The larger system includes, 
in addition to the carbons, the four (nonexchanging) hydrogens and the spin-1 nitrogen-14 
nucleus in alanine. The larger system's Hamiltonian has the form 

4 4 

Hs = Hs + nJ^^H ^H,z + ^ '^N ^N,z + ^Yl "^h/h ^h-^h + (22) 

4 3 3 4 

2 Zl I] -^c'h ^H,z^C,z + 2 5Z -^N.C ^N,z<,z + 2 5Z ^N,H ^N,z^H,z ' 

1=1 i=l i=l j = l 

where we have used the well-known fact that couplings between spins with distinct gyro- 
magnetic ratios can be truncated to just the secular (i.e. (TzCt^) part. 

The spin-1 nitrogen has an electric quadrupole moment and hence a short Ti, so its 
coupling to the other spins cannot be observed. As a result, the nitrogen can be omitted 
from the Hamiltonian, although it still plays a role in the relaxation of the other spins. The 
spin-lattice (Ti) relaxation time of the hydrogen atoms, on the other hand, is longer than 
the experiment, so here we must include the additional frequency shifts that depend upon 
their spin states. In other words, we take the (Th.z to be constants of the motion, and treat 
the carbons as an incoherent mixture of 2^ = 16 independent 3-spin systems, each with 
their resonance frequencies shifted by one of the 16 possible sums of ±Jc,h/2. Because we 
start with the spins in the high-temperature equilibrium state, each of these 16 independent 
evolutions contributes equally to the simulated density matrix. 

The computer search for the strongly modulating control sequences (pulses) is quite 
demanding, and it is important to keep the Hilbert space as small as possible during the 
associated simulations. For this reason the simulations described in this paper ignored the 
four hydrogen atoms in the alanine molecule, which were left aligned with the main magnetic 
field during the actual experiments but which nonetheless have couplings to the carbons on 

13 



the order of 150Hz. A better way, which we subsequently implemented and is described here 
for completeness, would be to average the results of 16 simulations over each of the 3-spin 
Hamiltonians 

^s-c = ^ E ('^c + E ^^ ^c;0^c,z + ^ E ^h/h ^h • ^'h , (23) 

for all combinations oi 6 = [Si, ■ . . ,5i\ = {0, 1}^, where the "c" stands for "control". This is 
much faster than diagonalizing the Hamiltonian with the four hydrogens included (16 times 
longer as opposed to ~ 16'^). 

During the free evolution periods between pulses, on the other hand, we applied additional 
(hard) 7r-pulses to the carbons to refocus this unwanted phase evolution, in accord with the 
free-evolution Hamiltonian used for these (much less demanding) simulations, 

4 3 

where the "f " now stands for "free" . The frequency evolution of the hydrogens can be ne- 
glected because they are left along the z-axis, thereby also avoiding population disturbances 
that would lead to nuclear Overhauser cross-relaxation and thereby memory effects. 

The search for strongly modulating control sequences was further simplified by assuming 
that the RF phase, amplitude, and frequency was piecewise constant, so the total unitary 
transformation was given by the product 

UL = n ^^ = n Texp(-i(tfc-tfc_i)if|, - i / dti/RF,fc(t)) (25) 

fc=l k=l -Jik-i 

Here, "T" denotes the usual time-ordering operator, and the external RF Hamiltonian during 
the /c-th interval is given by 

3 

Hn,At) = «,7ci?iE(^"^"*''"'^'^''^'^^c,xe'^"*''^'"^^^'^^-), (26) 

where at is the relative amplitude of the field during the fcth interval (cf. Eqn. [TT|) . This 
assumption allows the unitaries U^ for each interval to be calculated exactly by transforming 
to an interaction frame in which i^Rp, fc becomes time- independent, and diagonalizing the net 
Hamiltonian in that frame J3|. The spectrometer generates the control fields by applying a 
time-dependent voltage to a tuned resonator. The limitations on both the control circuitry 
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and the tuned resonator introduce time-dependent distortions of this modulation at the 
discontinuities between intervals. To account for this limitation, we monitor the field gener- 
ated in the control coil (the antenna that interacts with the spins) and pre-weight the time 
dependent wave-form to provide a close approximation to the desired shape. In addition, 
we use the measured modulation sequence in the simulator to follow the dependence of the 
propagator for small distortions. 

To achieve useable sensitivity, NMR is carried out on an ensemble of spatially distributed 
spins. The control field thus varies over the spatial extent of the sample, which is termed 
RF field inhomogeneity. This variation could be reduced by using a larger coil, a smaller 
sample, or by selecting a particularly homogeneous region of the sample via a magnetic field 
gradient. All of these options, however, reduce the signal-to-noise ratio. Fortunately, the 
spatial variation of the RF field is constant over time, and can be measured very accurately. 
This allows us to perform the simulations for each member of a histogram of the variations 
in RF amplitude across the sample, and to combine the results as an incoherent sum just 
as was described for the variations in the hydrogen spin states above. A total of 33 values 
was included in this RF inhomogeneity histogram (see Fig. 0}. 

The procedure used to find a modulation sequence that correctly implements any desired 
unitary operation seeks to maximize the fidelity between the desired unitary superoperator 
Uth ® ^th and the simulated superoperator Sop- The latter was obtained as an incoher- 
ent or Kraus operator over the 33 RF field strengths in the experimentally measured RF 
inhomogeneity histogram, i.e. 

33 33 

Sop(p) = 5^P.f/.pt/l = $^A£p4- (27) 

e=i 1=1 

Thus, the gate fidelity as defined in Eqn. ((7j) can be calculated directly from the Kraus 
operators A^ as 

Fe((f/,h ® f/th)%p) - wTm^r ! ^?r!frr\. = 2^ E tr((f/,h ® f/th)^^™ ® ^«.)) (28) 

^33 ,33 

= ^E^^((f^tU-)®(^tU.)) = ^Ei*^(^tUrn)i'. 

m=l m=l 

In simulations including the four hydrogen spins, the sum on the right-hand side of this 
formula must be increased by another factor of 16 in taking the partial trace over the 
hydrogen spins. 
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TABLE I: 

Measured carbon-13 decay rates for each type of product operator (s~^) 
in the alanine molecule used for the experiments. These include the 

inverses of the relaxation times Ti and T2 for all three carbon-13 spins. 
The "x|y" subscripts indicate that, in every operator of the corresponding 
product, all the labels to the left or to the right of the bar must be used 

together, thus indicating a total of two products with one expression. 



Type of Product Operator 



Decay Rate 



Type of Product Operator 



Decay Rate 



cr; 



0.032 ± 0.01 



cr; 



0.345 ± 0.01 



cr" 



0.583 ± 0.01 



CXzV.^ 



0.282 ± 0.01 



alaf 



0.458 ± 0.01 



aW. 



0.689 ± 0.01 



ala^a^ 



0.684 ± 0.01 



1 12 

'-''xly '''x|y'-'"z ' 
^x|y^z> f^xlyf^z^z 



1.89 ± 0.02 



2 12 

'-''xly '-'"z'-''x|y' 

<|yO"z' 0-z^x|yO-z 



3.19 ± 0.01 



_.3 _.1_3 

"x|y' "z"x|y' 

2 3 12 3 



1.68 ± 0.01 



12 123 

'^x|y^y|x' ^x|y^y|x^zi 



0"x|yC'"x 



|y' ^x|y^x|y^z 



6.93 ± 0.11 



13 123 

'-''x|y'^y|x' '-'"x|y'^z'^y|xi 

13 12 3 

'''xly'^xly '-'"x|y'^z'-''x|y 



3.56 ± 0.06 



2 3 12 3 

^\y y|x' z x|y y|x' 

^x|y^x|y ^z^x|y^x|y 



6.81 ± 0.10 



12 3 12 3 

'-^y|x'^x|y x|y '^x|y y|x'^x|y 

12 3 
'-^x|y'^x|y'-^y|x 



13.48 ± 0.81 



123 

"x|y xly^^xly 



14.58 ± 0.20 



The relaxation superoperator for the carbons in alanine was measured in the absence of 
RF fields, and found to have 98.5% of its norm along the diagonal in the product operator 
basis. This means that the various product operator components decay mono-exponentially, 
without cross-relaxation, so that it can be described by an 8 x 8 "Hadamard relaxation 



matrix" 



231 



3l|. This is shown pictorially in Fig. while precise values for the various 



types of rates seen in the figure are given in Table EJ 

To be complete, the relaxation superoperator should also be measured as a function 
of the applied RF fields (and include memory effects due to nuclear Overhauser effects 
with the hydrogens 32|). The QFT, however, efficiently mixes the states of the three 
carbon spins, so that on average over its implementation the decoherence is indistinguishable 
from a uniform, isotropic attenuation of all the product operator components. That is to 
say, QPT on the QFT is not able to provide any details regarding the physical relaxation 
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processes operative in alanine, save for an average overal rate of attenuation. Of course, 
QPT could be performed on a simpler gate, as in Ref. 15| where it was used to derive 
a relaxation superoperator for free precession, but the goal here is to use QPT to learn 
about the coherent and incoherent errors committed during implementation of a complex 
unitary transformation, the QFT. Therefore, this simple effective relaxation superoperator 
was assumed to accelerate the simulations of the overall QFT (no relaxation was assumed 
during the pulse design simulations, since they are not intended to correct for such effects). 

VI. AN OVERVIEW OF THE EXPERIMENTAL RESULTS 

Complete QPT was performed twice using different sets of strongly modulating RF control 
sequences to implement single spin rotations (Section IIV|) . The first such iteration was done 
with the sequences described in Ref. J3[, which refocussed the evolution under the alanine 
molecule's internal spin Hamiltonian [Hs in Eqn. ITn|l during the sequence. The second 
iteration used control sequences that not only refocussed the internal Hamiltonian, but also 
compensated for inhomogeneity in the RF field itself jlO|. Thus, the second set of control 
sequences was expected to produce much smaller incoherent errors than the first, while 
the coherent and decoherent errors were expected to be roughly the same for both. The 
improvement obtained with the RF-compensated set is illustrated in Fig. El which plots 
the correlations and attenuated correlations obtained for each of the 64 product operator 
basis states used as inputs for the two iterations against one another. Further evidence 
for substantial improvements with the RF-compensated control sequences may be obtained 
from the Kraus operator plots and statistics in Fig. [7| (see below for their interpretation). 
Having demonstrated this clear-cut improvement, all subsequent analysis will be given only 
for the results obtained wiht the compensated control pulses. 

The correlation provides an estimate of the accuracy of the experimental QFT implemen- 
tation, but without considering the loss of magnetization due to decoherence or incoherence. 
The average of the input state correlations over all 64 basis states was 0.96, while the aver- 
age correlation following application of the QFT was 0.82. A minimization search, however, 
found pure states which, after operating on them with the experimental superoperator, had 
a correlation with the same pure state following the theoretical QFT as low as 0.45. The 
average attenuated correlation (or gate fidelity 'j^), which also takes into account the loss of 
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magnetization, was 0.64, indicating that about 22% of the magnetization was lost over the 
ca. 30 ms needed to implement the QFT (see Section HVJ) . These numbers are in-line with 
expectations based on other recent apphcations of the RF-compensated control sequences 
j^. About half of this magnetization loss was expected due to intrinsic decoherence, and 
since approximately the same amount of decoherence occurred during the input states' read- 
out as during the output states', the remainder is probably due to residual uncompensated 
incoherence and/or imperfect decoupling of the protons from the carbons used as qubits. 
These issues will be discussed in more detail in later sections. A complete list of all the 
output state (attenuated) correlations may be found in Fig. |H] below. 

Rather than looking at the action of the QFT on states, one can also look directly at 
the Kraus operator sum computed from the completely positive part of the experimental 
supermatrix Mobs (Section IIII|) . Each Kraus operator Ak has an associated amplitude, 
o-k = ||^fe||/"\/8 = A/Ajt/8 (where A^ is the kth eigenvalue of the associated Choi supermatrix), 
and since in a perfect implementation the desired unitary operator f/qpT would be the only 
Kraus operator with nonzero amplitude, one expects the Kraus operator with the largest 
amplitude to be at least fairly similar to t/QPT- This is confirmed by the plots of the real 
part of the largest Kraus operator, shown in Fig. [71 which had an amplitude ai = 0.86 
anda correlation with the real part f/gpT of 0.95, implying a net coherent error of roughly 
5%. The second largest Kraus operator, also shown, had an amplitude of 02 = 0.34, and 
was also rather close to unitary although uncorrected with J/qft- We expect it to be a 
rough approximation to the largest unitary operator in any sum of unitary transformations 
making up the incoherent error. Finally, there are another 32 essentially non-unitary Kraus 
operators in the completely positive part of the experimental QFT supermatrix with smaller 
amplitudes (Fig. [7j). 

The non-completely positive part of the QFT supermatrix, which is obtained from the 
eigenvectors associated with the negative eigenvalues of its Choi supermatrix, will also be 
of interest in what follows. The positivity of the QFT supermatrix, as defined at the end 
of Section ITTTl was only 0.60, but the ratio of the smallest to the largest Choi supermatrix 
eigenvalues was —0.075, indicating that the negative eigenvalues were rather small in mag- 
nitude in comparison to the positive. The ratio with the second largest was —0.48, and 
the third had almost the same magnitude as the smallest. These observations, together 
with the well-known sensitivity of the eigenvectors of nearly degenerate eigenvalues to small 
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perturbations in the elements of the matrix from which they come, imply that rather little 
information about the errors made in the QFT implementation can be gleaned from the 
individual Kraus operators after the first two. Only the subspaces spanned by the eigenvec- 
tors associated with all of the smaller positive, or perhaps negative, eigenvalues is likely to 
be statistically significant. 

VII. COMPARISON WITH THE MODEL 

To assess the precision and completeness of our system-plus-apparatus model, the com- 
plete set of experiments involved in QPT on the QFT was simulated using the mathematical 
model of the system-plus-apparatus described in Section El The compatibility of the simu- 
lated results with the experimental provides a rigorous test for the accuracy of the model, 
which takes all the known significant imperfections of the experimental apparatus into ac- 
count. All of the unitary operations performed during these simulations were implemented 
using exactly the same strongly modulating control sequences that were used for the exper- 
iments, and complete state tomography was performed for each input and output state. It 
should be noted, however, that these states were reconstructed directly from the observables 
in the simulated density matrices, without further simulating the spectra and fitting them 
as required in experimental state tomography (which of course gives rise to additional errors 
in the actual experiments). Figure |H1 shows the correlation and attenuated correlation of the 
initial and final states obtained from these simulations. 

Figure M plots the sorted Kraus operator amplitudes obtained from the simulated su- 
permatrix, along with those from the experimental supermatrix for comparison (cf. Fig. [7j). 
The negative values plotted are actually the negative square roots of the corresponding Choi 
matrix eigenvalues, and are shown to illustrate that the experimental supermatrix was sig- 
nificantly further from being completely positive than was the simulated (the positivity of 
the simulated was 0.86, as opposed to 0.60 for the experimental). The most likely reason 
for this is the absence in the simulations of the additional errors expected from fitting the 
spectra to extract the product operator amplitudes for state tomography. Also shown once 
again is the real part of the theoretical QFT unitary matrix (cf. Fig. Ej), along with the real 
parts of the matrices of the largest Kraus operators from the experimental and simulated 
supermatrices. Finally, the corresponding best unitary approximations to the largest Kraus 
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operators, obtained by setting their singular values to unity, have their real parts shown on 
the bottom right-hand side of Fig. IHl It may be seen that there is a good correspondence 
between the largest Kraus operators as well as between their best unitary approximations, 
with a correlation between the simulation and the experiment of 0.90 in both cases. As noted 
previously, the smaller Kraus operators cannot be expected to correspond significantly to 
one another. 

Since the Fourier transformation converts between the position and momentum bases, 
it has a very simple interpretation in phase space. Classically, a Fourier transform rotates 
phase space by 90°. In quantum phase space, the QFT superoperator has only one nonzero 
element equal to unity in each row and column, and so constitutes a permutation matrix 
J33[. This is shown in Fig. ^[ along with the corresponding plots for the simulated and 
experimental supermatrices. One can see immediately from these plots that there are errors 
in the QFT implementation. The phase space basis, however, consists of operators that 
are neither Hermitian nor products of any underlying Hilbert space basis, and which are 
not related in any simple way to the physical operators used to implement the QFT. Thus, 
although the action of the QFT is very easy to understand in the phase space basis, it is 
very hard to interpret the discrepancies between the theoretical superoperator and observed 
supermatrix in terms of implementation errors in this basis. It is worth noting, nonetheless, 
that in the phase space basis the errors appear more white than Gaussian especially in the 
simulation. 

Other bases, in which the discrepancies between the simulated and experimental super- 
matrices are also clearly manifested, include the computational basis (or Zeeman basis; see 
Fig, nil) , and the product operator basis (shown in Fig. IT^ . where the latter consists of all 
possible products of the Paul, ,„at*es ... ., and .. of different spins R (see Section O 
above). It may be observed that the theoretical superoperator contains several fixed points 
in the product operator basis, which provide us with another interesting metric for the pre- 
cision of our implementation. Two of these fixed points are cr^cr^ and {a^ + cr^)/2. The 
correlations between these fixed points, before and after applying the experimental and sim- 
ulated superoperators, are 0.91 and 0.93 respectively. Further information can be obtained 
from plots of the individual rows, which depict how much each input state contributes to a 
given output state, as shown in the relatively simple case of the cr^cr^ fixed point in Fig. [T3I 

The greatest part of the deviation between the simulated and experimental results is 
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expected to be due to the propagation of small coherent errors. Thus it is interesting to 
compare the simulated and experimental supermatrices after correcting for these errors as 
completely as possible. This is most simply done by left-multiplying by the inverse of the 
superoperator obtained from the best unitary approximation to the largest Kraus operator, 
Ui ® Ul, followed by the theoretical unitary superoperator J/qft ® Uqpj, i.e. 

3VCcor,obs|sim = (f^QFT ® f/gFj) (f^l,obs|sim ® f^l,obs|sim) Mobs|sim- (29) 

The resulting supermatrices in the phase space basis are displayed in Fig. El along with 
the exact QFT supermatrix in the same basis for comparison. It may immediately be seen 
that there has been a considerable improvement in the similarity of these matrices, as is fur- 
ther confirmed by a correlation between the corrected simulated and corrected experimental 
supermatrices of 0.94, between the corrected simulated and theoretical of 0.99, and between 
the corrected experimental and theoretical of 0.95. The correlations between the simulated 
and experimental correction factors, i.e. the products of the unitary superoperators in the 
above equation, are however only 0.90, indicating that the cumulative effects of coherent er- 
rors that were not taken into account by the simulations over the course of the experiments 
was roughly 10%. It should further be noted that whereas the experimental supermatrix is 
significantly more strongly correlated with the theoretical than it is with the simulated, after 
they have been corrected their correlations and attenuated correlations with the theoretical 
are very nearly the same. 

VIII. DISCUSSION 

We have shown that our model of the system-plus-apparatus is able to predict many 
details of the experimental results (see Table |n] for all the correlations among the theoret- 
ical, simulated and experimental supermatrices, together with the attenuated correlations 
to the theoretical). Specifically, simulations based on the Hadamard relaxation operator 
shown in Fig. and Table H] have allowed us to establish that the results contain no specific 
information on the decoherence rates and processes operative in our system, since these are 
averaged by the complex sequence of transformations that make up the QFT. Their net 
effect can therefore be modeled as a simple uniform attenuation of all the product operators 
in the density matrix other than the identity. This is explicitly demonstrated by the close 

21 



TABLE II: 

The correlation coefficients between all pairs of super matrices, labelled as described in 

the main text, together with the attenuated correlations (as defined in section II) 

between each and the theoretically exact supermatrix of the quantum Fourier transform. 



Supermatrix 


Correlations Among All Pairs of Supermatrices 


Attenuated 
Correlation 


Theoretical 


1.00 0.89 0.99 0.82 0.95 


1.00 


Simulated 


1.00 0.90 0.79 0.85 


0.68 


Corrected Sim. 


1.00 0.81 0.94 


0.76 


Experimental 


- - - 1.00 0.87 


0.64 


Corrected Exp. 


- - - - 1.00 


0.74 



correspondence between the supermatrix eigenvalues shown in Fig. ^J both from simula- 
tions using the Hadamard relaxation operator (red " * " ) as well as from simulations taking 
no account of relaxation save by scaling down the nonidentity components of the density op- 
erators by a factor of 0.82 (blue "o"). In general, of course, relaxation cannot be accounted 
for by a single attenuation factor, and in fact methods similar to those described here have 
jeen used to determine the complete NMR relaxation superoperator of a two-spin system 
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Figure ITHl also shows the eigenvalues of the supermatrix obtained from simulations which 
included incoherent errors from RF field inhomogeneity, first without taking relaxation or 
readout errors into account (cyan "•"), and second from simulations which included input 
and output state readout and took relaxation into account by scaling the nonidentity compo- 
nents of the density matrices read out by 0.82 (green "x"). These plots show rather clearly 
that most of the additional dispersion seen in the simulations including incoherent errors 
stems directly from those errors, and was not an unintended consequence of the additional 
errors made during the readout operations needed for QPT (recall that incoherent errors 
were included in simulating the readout operations). 

Nevertheless, the spread in the eigenvalues of the experimental superoperator (see Fig. El 
below) were significantly larger than those in the simulations of the full QPT procedure, 
and it was not possible to match them in a one-to-one fashion. In part this may be due 
to systematic errors in the least-squares fitting procedure by which the observables (nixx 
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in Fig. 12)) were extracted from the NMR spectra, in particular errors in distinguishing the 
absorptive and dispersive peak components. Such errors are difficult to simulate accurately 
and therefore not included in our simulations of the full QPT procedure, but are consistent 
with the fact that the largest single-spin corrections were z-rotations (as described below). 
Another, probably more important reason for these apparent discrepancies between the 
model and the experiments lies in the propagation of many small coherent errors over the 
course of the QFT implementation. 

Even though the simulated single-qubit gate fidelities were all better than 0.99, they 
were also less than unity. The main thing that limited these fidelities was that the number 
of parameters, and hence the number of time intervals within which the RF amplitude 
and phase was held constant, had to be kept as low as possible during the optimizations by 
which the strongly modulating pulses were designed. Although the cumulative effects of such 
small errors could become significant, since they are included in the simulations they are not 
likely to be the source of the additional eigenvalue dispersion that is seen in the experimental 
results. A more important source of coherent errors is expected to be due to the fact that the 
protons were not included in any of the simulations, and their couplings to the carbons used 
as qubits will show up as phase rotations. In addition, the model ignores such fine details of 
the apparatus as the limits on the rise-and-fall times of the transmitter used to generate the 
RF control sequences, and more generally its frequency response characteristics. Although 
the simulation of reactive circuits is more complex than the simulation of resistive circuits, 
in due course we will also include such effects in the model, and continue to refine it until 
the only remaining discrepancies between the simulated and experimental results lie in the 
intrinsic measurement errors. 

Unwanted bilinear interactions are refocussed during the QFT implementation by the 
application of vr-rotations to the individual spins, using traditional NMR methods |32j |. 
Errors in these rotations will contribute to the bilinear errors only to second order, so these 
bilinear errors depend mainly upon the timing of the pulses, which can be controlled very 
accurately. As a consequence, we expect that the residual coherent errors that were not 
taken into account by the simulations will be primarily single-spin rotations. To test this 
hypothesis, the error operator Ua = f^i,exp^isim ^^-s taken, where t/i,exp is the best unitary 
approximation to the largest Kraus operator Ai of the experimental supermatrix (cf. Fig. Ej), 
and f/i,sim is similarly the best unitary approximation to the largest Kraus operator of the 
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simulated superinatrix. Using a numerical search, the product of three single-spin rotations 
was found that fit Ua best in the least squares sense. The resulting unitary U^ ® W^ ® U^ 
had a correlation coefficient with Ua of 0.96, in support of the hypothesis. 

The angles and directional cosines of the axes of these three rotations are shown in Table 
11111 while the eigenvalues of the resulting corrected superoperators are shown in Fig. ^1 
It should be noted that, although the correlation between Ui^exp and f/i,sim increased from 
0.90 to 0.96 on left-multiplying Ui by the Hermitian conjugate of this product of single-spin 
rotation operators, these rotations are the cumulative result of many small rotation errors 
and are not simply traced back to any single short-coming in the experiments. For com- 
pleteness, the axes and angles of the single-spin rotations that best fit the error operator 
Ua = Ul 5jj^ f/i,exp are also shown in Table 11111 together with the eigenvalues of the corre- 
sponding corrected superoperators in Fig. ^| The product of these rotations similarly has 
a 0.97 correlation with Ua, but in this case one must right-multiply Ui^exp by the Hermitian 
conjugate of the product to correct it. The close co-incidence between the angles of rotation 
about the x-axis on spin 1 in the first case and about the z-axis on spin 3 in the latter case 
is expected, since ^qftcTx^qft — ^l (recall cr^o"^ is a fixed point of the QFT). 

Finally, it is of interest to demonstrate that despite a substantial number of negative 
eigenvalues in the Choi matrix of the experimental supermatrix, it is not necessary to change 

TABLE III: 

The two triples of single-spin rotations that optimize the correlation coefficient 

between the best unitary approximations to the largest Kraus operators of the 

experimental and simulated supermatrices, when the experimental supermatrix is left 

or right multiplied by the unitary matrix corresponding to each one of these triples. 



Side 


Spin 


X, y &: z-Directional Cosines 


of Rotation Axis 


Rotation Angle 


left 


1 


-0.992 0.007 


0.123 


36.7° 


left 


2 


-0.386 -0.243 


0.890 


9.0° 


left 


3 


0.703 0.701 


-0.123 


16.9° 


right 


1 


0.059 -0.031 


0.998 


14.2° 


right 


2 


0.227 -0.512 


0.829 


10.2° 


right 


3 


-0.092 -0.292 


-0.952 


38.3° 
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it much in order to obtain a superniatrix which represents a completely positive and trace- 
preserving superoperator. For this reason the superniatrix Mcptp which best-fit the Choi 
matrix of the experimental superniatrix subject to the constraint that it was both positive 
semidefinite and satisfied the trace-preservation conditions was computed as described in 
Section 11111 Although this procedure made essentially no change in the largest Kraus op- 
erator (as expected), it did have a significant effect on the experimental superniatrix as a 
whole. The correlations between this CPTP-fit and the other supermatrices that we have 
dealt with up to now are given in Table IIVI along with those to the original experimental 
superniatrix for comparison. This shows that even though imposing the complete positivity 
constraint on the experimental observations did not change the superniatrix very much, the 
change was distinctly in the right direction since it improved the correlation with both the 
simulated and theoretical supermatrices. This is further confirmed by Fig. El which shows 
the eigenvalues of Mcptp, along with those of the superoperator f7i,cPTP ® f^i,CPTP obtained 
from the best unitary approximation to its largest Kraus operator (scaled down so as to have 
the same trace as Mcptp), and those of the experimental supermatrix Mobs for comparison. 
From this we see that the eigenvalues of Mcptp are closer to being cocircular, indicating that 
it is closer to an attenuated unitary than was Mobs, but that only after taking the unitary 
part did they become perfectly cocircular. 



TABLE IV: 

The correlation coefficients between Theoretical, Simulated and Experimental supermatrices 
(see text) and the optimum completely positive and trace-preserving approximation Mcptp to 
the experimentally determined supermatrix Mobs, as well as the best unitary approximation 
to its largest Kraus operator C/i^cptp "X" f7i,CPTP and Mobs itself for comparison. 



Supermatrix 


Correlation with 

f^l,CPTP ® f/l,CPTP 


Correlation with 
Mcptp 


Correlation with 

Mobs 


Theoretical 

Simulated 

Experimental 


0.86 
0.82 
0.95 


0.89 
0.82 
0.97 


0.82 
0.79 
1.00 
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IX. CONCLUSIONS 

In conclusion, we have implemented quantum process tomography of the quantum Fourier 
transform on a three-qubit NMR quantum information processor. The overall gate fidelity 
(attenuated correlation between superoperator matrices) was 0.64, whereas the (unattenu- 
ated) supermatrix correlation was 0.82 (see Table HT]) . Judging by the fact that making the 
unitary part of the largest Kraus operator correspond as closely as possible to the theoreti- 
cal QFT gave a correlation of 0.95, we conclude that the loss of fidelity due to incoherence 
and/or measurement errors during state tomography was of order 5%. The loss of mag- 
netization due to incoherence and decoherence, on the other hand, was 0.64/0.82 = 0.78, 
i.e. about 22%. This implies that the cummulative effects of coherent errors reduced the 
fidelity by about 0.82/0.95 = 0.86, or 14%, consistent with the fact that the correlation 



between the largest Kraus operator and the theoretical QFT unitary was 0.925 ~ \nL86. 

More importantly, QPT of the QFT has enabled us to validate the essential correctness 
of our model of the system-plus-apparatus used for the experiments in great detail, and 
to isolate its remaining shortcomings. It has further prompted us to develop a range of 
data analysis and visualization techniques for quantum process tomography, which should 
be broadly applicable in quantum information processing. The experiments described here 
demonstrate the precision with which complex quantum dynamics can be controlled, and 
highlights the significance of liquid-state NMR as a test bed for achieving such control. While 
full QPT on larger quantum systems will never be practical, the analysis done here should 
serve as an initial guide as to how information about the errors in quantum information 
processors can be extracted, and perhaps someday, how to debug a quantum computer. 



26 



ACKNOWLEDGEMENTS 

The authors thank J. P. Paz for helpful discussions. This work was supported by 
ARDA/ARO grants DAAD19-01-1-0519 & DAAD19-01- 1-0654, DARPA grant MDA972- 
01-1-0003, NSF grant EEC 0085557, the Air Force Office of Sponsored Research, and by the 
Cambridge-MIT Institute, Ltd. 



[1] B. Schumacher, Phys. Rev. A, 54, 2614, (1996). 

[2] G. Teklemariam, E. M. Fortunate, M. A. Pravia, T. F. Havel, and D. G. Cory, Phys. Rev. 

Lett. 86, 5845 (2001). 
[3] E. M. Fortunato, M. A. Pravia, N. Boulant, G. Teklemariam, T. F. Havel, D. G. Cory J. 

Chem. Phys. 116, 7599-7606, (2002). 
[4] N. Boulant, K. Edmonds, J. Yang, M. A. Pravia, and D. G. Cory Phys. Rev. A 68, 032305 

(2003). 
[5] 1. L. Chuang, M. A. Nielsen, J. Mod. Optics 44, 2455 (1997). 
[6] J. F. Poyatos, J. I.Cirac, P. Zoller, Phys. Rev. Lett. 78, 390 (1997). 

[7] G. M. D'Ariano, L. Maccone, Phys. Rev. Lett. 80, 5465 (1998); G. M. D'Ariano, M. G. A. 
Paris, M. F. Sacchi, Advances in Imaging and Electron Physics (Academic Press, Englewood 
Cliffs, NJ), in press, http:arxiv.org/abs/quant-ph/0302028. 
[8] Y. S. Weinstein, M. A. Pravia, E. M. Fortunato, S. Lloyd, D. G. Cory, Phys. Rev. Lett. 86, 

1889 (2001). 
[9] M. H. Levitt, R. Freeman, J. Magn. Reson. 33, 473 (1979); M. Levitt, J. Magn. Reson. 48, 234 
(1982); R. Tycko, Phys. Rev. Lett. 51, 775 (1983); A. Shaka, R. Freeman, J. Magn. Reson., 
55, 487 (1983); M. Levitt, Prog. Nucl. Magn. Reson. Spect., 18, 61 (1986). 
[10] M. A. Pravia, N. Boulant, J. Emerson, A. Farid, E. M. Fortunato, T. F. Havel, R. Martinez, 

D. G. Cory J. Chem. Phys. 119, 9003 (2003). 
[11] K. Kraus, States, Effects, and Operations (Springer- Verlag, Berlin, FRG, 1983). 
[12] D. W. Leung, J. Math. Phys. 44, 528 (2003). 
[13] T. F. Havel, J. Math. Phys. 44, 534 (2003). 

[14] A. M. Childs, I. L. Chuang, D. W. Leung, Phys. Rev. A 64, 012314, (2001). 
[15] N. Boulant, T. F. Havel, M. A. Pravia, D. G. Cory Phys. Rev. A 67, 042322 (2003). 

27 



[16 

[17: 

[18 

[19 
[20 
[21 

[22; 

[23 
[24 
[25; 
[26 

[27; 

[28 



P. W. Shor, Siam. J. Comput. 26 1484-1509 (1997). 
R. Jozsa, Proc. Roy. Soc. Lond. 454 (1998). 



C. Zalka, Proc. Roy. Soc. Lond. A 454 313-322, (1998); S. Wiesner, |^iiaiit^h/9603028. 
R. Schack, Phys. Rev. A 57 (1998). 

D. Coppersmith, IBM Research Report RC19642 (1994). 
M. D. Price, T. F. Havel and D. G. Cory, New J. Phys. 2, 10 (2000). 
Y. S. Weinstein, S. Lloyd, J. Emerson, D. G. Cory, Phys. Rev. Lett. 89, 157902 (2002). 
T.F. Havel, Quant. Inf. Proc, 1, 511, (2003). 
P. Pechukas, Phys. Rev. Lett. 73, 8, (1994). 
P. Stelmachovic, B. Buzek, Phys. Rev. A 64, 062106, (2001). 
P. Stelmachovic, B. Buzek, Phys. Rev. A 67, 029902, (2003). 
N. Boulant, S. Furuta, J. Emerson, T. F. Havel, D. G. Cory, l quant-ph/0312116[ 
D. G. Cory, R. Laflamme, E. Knih, L. Viola, T. F. Havel, N. Boulant, G. Boutis, E. Fortunato, 
S. Lloyd, R. Martinez, C. Negrevergne, M. Pravia, Y. Sharf, G. Teklemarian, Y. S. Weinstein, 
Z. H. Zurek, Fortsch. Phys. 48, 875 (2000). 

[29] J. A. Jones, Fortsch. Phys. 48, 325 (2001). 

[30] S. S. Somaroo, D. G. Cory, T. F. Havel, Phys. Lett. A, 240, (1998). 

[31] T.F. Havel, Y. Sharf, L. Viola, D.G. Cory, Phys. Lett. A, 280, 282 (2001). 

[32] R. R. Ernst, G. Bodenhausen, A. Wokaun, Principles of Nuclear Magnetic Resonance in One 

and Two Dimensions, Oxford Univ. Press, Oxford (1994). 
[33] C. Miquel, J.P. Paz, M. Saraceno, Phys. Rev. A, 65, 062309, (2002). 



28 



FIGURE CAPTIONS 

FIG. 1: Circuit diagram for implementation of the quantum Fourier transform for three 
qubits. For each qubit j, starting with the most significant, a series of conditional phase 
gates are implemented between qubit j and all qubits more significant than j, followed by a 
Hadamard (H) on j. The amount of phase added is 9jk = n/2^^''. A bit reversal (two-headed 
arrow) completes the QFT. 

FIG. 2: A schematic of the QPT implementation used in this paper, together with two 
equations describing how the readout operations Uro and measurements mxx (xx = in, op) 
are obtained from the ensemble of quantum systems described by the density matrix p 
(cr_ = |((Tx — io'y))- A set {Sin} of (not necessarily unitary) operations is performed, each in 
a separate experiment, on the equilibrium state peq of the spin system to create a complete 
set {pin} of input states. Each input state is then determined by quantum state tomography, 
i.e. by repetition of the experiment with different readout pulses, {IXro}, appended to each 
repetition. The readout pulses rotate unobservable components of the density matrix into 
observable components, the mean values {mjn} of which are equal to the unobservables' mean 
values before rotation. Measurement of these mean values thus allows for the reconstruction 
of the input state that was actually created. Next, the set of input states is recreated one 
state at a time, and the operation §op of interest (in this paper the QFT) applied to each 
in turn. This gives a complete set of output states, {pop}- Once again, readout of each pop 
requires that the experiment be repeated followed by different readout operations {IXro}- The 
mean values of the complete set of observables {mop} that is measured allows reconstruction 
of the set of output states actually created. Finally, these estimates of {pop}, together with 
the earlier estimates of {pin}, are used to reconstruct an estimate of the operation §op via 
Eqn. 0. This reconstructed operation is not expected to be unitary due to decoherence 
during the application of §op as well as errors in the state tomography procedure. Also, 
due to incoherent errors across the ensemble involved in NMR experiments, the supermatrix 
Mobs estimated via this QPT procedure may not correspond to a completely positive, trace 
preserving superoperator. 

FIG. 3: Eigenvalues of simulated supermatrices in the complex plane. The main figure 
shows the location of the eigenvalues with respect to the unit circle, while the subplot shows 
the eigenvalue locations in greater detail. This simulation was done for alanine (see text). 
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and includes both coherent errors due to imperfect pulses and incoherent errors due to the 
inhomogeneous RF power across the sample. The eigenvalues indicated by circles (o) are of 

1 9 1 1 

the product Sigox Sgoy of ^^e supermatrix describing an evolution SgQy of a 90° pulse on spin 
1 about the negative y-axis times a second supermatrix Sigox for a 180° rotation of spins 1 
and 2 about the x-axis. The eigenvalues indicated by crosses (x) are for the supermatrix 
Sboth describing the net evolution after concatenating these pulses together in the simulator. 
Due to the presence of incoherence the eigenvalues are not the same. In particular, the 
presence of incoherence leads to an apparent decoherence, as evidenced by a reduction in 
the eigenvalue magnitudes. The average eigenvalue reduction is 1.2% for Sboth and 0.9% for 

c?l,2 Pi 
"I80x*^90y- 

FIG. 4: The 33 point RF profile shows how much of the sample sees what fraction of the 
desired power (so 1.0 on the horizontal axis is the desired power). This profile was measured 
experimentally and used for the simulations. 

FIG. 5: Measured Hadamard relaxation operator versus the product operator basis for the 
three carbons of alanine in the absence of RF fields, with rates color-coded as indicated by 
the legend. The labels on the axes are a short-hand for the product operators associated 
with with X and y directions on the three carbons, with e.g. XIX ^^ cr^cr^- The labels 
for each entry of the matrix are obtained by analogy with the product rules for the Pauli 
operators; for example, the matrix entry in the row labeled IXX and column labeled YIY is 
YXZ ^^ a^a^a^ 123:] • See Table HI for the numerical values of the various entries. 
FIG. 6: The various correlations (see Section |n]) for the 64 states of the product operator 
basis before and after applying the QFT, where those obtained from the uncompensated 
strongly modulating control sequences have been plotte d ag ainst those obtained with the 
sequences designed to compensate for RF inhomogeneity [10|. To spread the points out, the 
(co)ordinates plotted are log^ol^ + (1 ~^))5 although the axes have been labeled by the actual 
correlations < a; < 1. Circles (o) are used for the initial state correlations, while triangles 
(a) denote the final state correlations and stars (*) indicate the final state attenuated 
correlations. For the uncompensated sequences the average initial state correlation is 0.93, 
the average final state correlation is 0.64, and the average attenuated correlation is 0.37. 
For the compensated sequences, the average initial state correlation is 0.96, the average 
final state correlation is 0.82, and the gate fidelity 0.64. Moreover, since a large majority 
of the points lie below the diagonal in the plot (dashed line), the correlations with the 
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compensated pulses were in fact better in almost every case. 

FIG. 7: The bar graph on the left shows the amplitudes ak = \\Ak\\/V8 of experimental 
QFT Kraus operators, where the light grey portions of the bars are the amplitudes obtained 
from QPT done with the uncompensated control sequences, and the dark grey bars those 
obtained from QPT done with the RF-compensated control sequences (see text). With the 
uncompensated control sequences the dominant Kraus operator had an amplitude a" = 0.72, 
whereas the compensated had a\ = 0.86. The second largest Kraus operators had amplitudes 
of a2 = 0.30 and 03 = 0.34, respectively, after which the amplitudes decreased much more 
slowly. The plots on the right (which are shaded by column number to aid viewing) compare 
the real part of the desired unitary operator J/qft (top center) to the real part of the largest 
Kraus operator Ai obtained using the uncompensated and compensated control sequences 
(left and right-hand plots on the second row, respectively), and to the corresponding second 
largest Kraus operators A2 (bottom row). The correlations between real parts of these 
operators were C{^UQfT,^A"^) = 0.78, C{^UQf:j,^Al) = 0.95, C(3f?f/QFT, ^^^2) = 0.06 and 
C(3fJf/QFT, 3f?A^) =0.02. 

FIG. 8: The differences between the (attenuated) correlations of the output states and 
their mean values are displayed for each input state, which are indicated by the symbols on 
the left vertical axis with e.g. IXY ^^ o'x^y- The simulated-vs-theoretical state correlations 
are in blue, and the experimental-vs-theoretical state correlations are in yellow, while the 
corresponding attenuated correlations are in cyan and brown, respectively. The correspond- 
ing mean values of the (attenuated) correlations were 0.89, 0.82, 0.83 and 0.63, respectively. 
The pair of dashed vertical lines about the center mark one standard deviation for the input 
state correlations (not shown), which are expected to be determined primarily by measure- 
ment errors so that any values within these bounds are certainly almost entirely due to 
noise. The larger deviations in the correlations are due to the propagation of coherent error 
(see text), while the yet larger deviations in the attenuated correlations include the effects 
of decoherence (modeled as a uniform attenuation in the simulations). 
FIG. 9: The left-hand plot shows the real amplitudes, of experimental QFT Kraus oper- 
ators (dark grey) and those obtained from the QPT simulation based on the system-plus- 
apparatus model (light grey), where the negative values plotted are actually the negative 
square roots of one-eighth the absolute values of the corresponding Choi matrix eigenvalues. 
It may be seen that the experimental deviated significantly more from being completely 
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positive than the simulated, most hkely to the absence of errors from fitting the spectra 
in the latter (see text). The dominant Kraus operator has an amplitude 0.86 for the ex- 
periment and 0.87 for the model. The right side of the figure compares the real part of 
the theoretical unitary QFT operator (top-middle) to the largest Kraus operator from the 
experiment (middle-left) and simulated (middle-right) supermatrices. Also shown are the 
real parts of the corresponding best unitary approximations to these largest Kraus operators 
(bottom line). The correlation between the two largest Kraus operators, and between the 
corresponding best unitary approximations, was 0.90 in both cases. 

FIG. 10: The theoretical (left), simulated (middle) and experimental (right) supermatrices 
of the QFT in the phase space basis, wherein all their elements are necessarily real. Since 
the theoretical is simply a permutation matrix in this basis, the differences between it and 
the other two are easily seen, but the physical meanings of these differences are obscure. 
For example, the level of "background noise" is noticeably larger in the simulation than in 
the experiment, although both have the maximum entry of each row (column) in the same 
place as the theoretical. 

FIG. 11: The real parts of the theoretical (left), simulated (middle) and experimental 
(right) supermatrices in the computational (Zeeman) basis, where the theoretical superop- 
erator is given in terms of the unitary operator by U^^^ U^ . The correlations and 
attenuated correlations among these matrices may be found in Table |n] Although the pat- 
tern of the theoretical superoperator's matrix elements is visible in both the simulated and 
the experimental supermatrix, it is still not at all obvious from this vantage point what 
errors have occurred during implementation or how they can be corrected. 
FIG. 12: The theoretical (left), simulated (middle) and experimental supermatrices in the 
product operator basis, wherein all the elements are necessarily real (see text; the scale has 
been reduced by a factor of 2 to create more contrast). The correspondence between the 
rows / columns of this supermatrix and the individual product operators is the same as that 
given by the labels on the vertical axis of the bar graph in Fig. |H1 Because all the entries 
of the corresponding "real" density matrix are the expectation values of observables which 
transform nicely under rotations 23|, the physical interpretation is somewhat easier in this 
basis. This is manifest, for example, in the appearance of several fixed points with a 1 on 
the diagonal and zeros elsewhere in the same row (see text). 
FIG. 13: Theoretical (left), simulated (middle) and experimental (right) supermatrices in 
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the phase space basis after correcting the latter two by a unitary supermatrix designed to 
bring them as close as possible to the theoretical (as described in the text). The levels of 
background noise are a great deal less than in the corresponding uncorrected supermatrices 
(Fig. CUD. 

FIG. 14: As may be seen from row 14 of the supermatrices shown in Fig. ^1 the operator 
crlcrl is a fixed point of the QFT. The upper bar graph displays the corresponding row of 
the theoretical (dark grey), simulated (medium grey) and simulated with coherent errors 
corrected (light grey) supermatrix elements versus the product operator basis, whereas the 
lower bar graph displays these same statistics for the experimental supermatrix. In contrast 
to the phase space and Zeeman basis (not shown), there is a significant correlation between 
the experimental and simulated values; specifically, the simulated-to-experimental correla- 
tion coefficients were 0.80 and 0.94 before and after correction, respectively. In the present 
case these correlations are due almost entirely to the single large diagonal value, but similar 
correlations were obtained for all other pairs of rows; the average correlations were 0.78 and 
0.94 before and after correction, respectively. 

FIG. 15: Plot of the eigenvalues of simulated QFT superoperators in the complex plane, 
specifically, the results of simulations with only coherent errors (blue "o"), the results in- 
cluding relaxation using the relaxation superoperator shown in Fig. Eland Table|l](red "* "), 
the results including incoherent errors but no relaxation (cyan "•"), and the same following 
simulation of the full QPT procedure (green "x"), wherein decoherence was mimicked sim- 
ply by scaling down the simulated QFT superoperators following the state readout pulses 
by a factor of 0.82 (save for the component with eigenvalue 1, which is needed for trace- 
preservation). To aid the comparison, the eigenvalues of the QFT superoperator including 
only the coherent errors (blue "o") were scaled to have the same RMS value as those ob- 
tained with the simulated relaxation, thereby showing that the net effect of relaxation on 
the results of the QFT is simply to lower the amplitude of the eigenvalues without otherwise 
dispersing them in the complex plane. It may also be seen that the effect of incoherent 
errors on the eigenvalues is both to increase their angular spread and scale them down pro- 
portionately, so that they tend to move along arcs inside of and tangent to the unit circle, 
and that the additional errors introduced by the tomography procedure do not alter them 
greatly. 
FIG. 16: The blue squares (n) are the eigenvalues of the simulated supermatrix Msim, 
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including coherent, incoherent and readout errors, while the magenta triangles (a) are those 
after correcting Mobs by post-multiplication with the product of single-spin rotations U^ ® 
f/^ ® U^ and the orange diamonds after pre- multiplication with a product of different single- 
spin rotations (see Table ITTT|l . Although we cannot quite unambiguously match up pairs of 
eigenvalues, it is clear that they have been made very similar by these corrections. 
FIG. 17: The blue plus-signs (+) are the eigenvalues of the experimental supermatrix Mobs, 
while the red five-poined stars (*) are those its best completely positive and trace-preserving 
approximation Mcptp, and the green six-pointed stars (♦) are those of the superoperator 
corresponding to the best unitary approximation to its largest Kraus operator t7i,cPTP ® 
^i,CPTP (see text), scaled down to have the same trace as Mcptp- 
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